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We explore the stability properties of multi-field solutions of assisted inflation type, where several 
fields collectively evolve to the same configuration. In the case of noninteracting fields, we show 
that the condition for such solutions to be stable is less restrictive than that required for tracking 
in quintessence models. Our results, which do not rely on the slow-roll approximation, further 
indicate that to linear order in homogeneous perturbations the fields are in fact unaware of each 
other's existence. We end by generalizing our results to some cases of interacting fields and to other 
background solutions and dynamics, including the high-energy braneworld. 
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I. INTRODUCTION 

The dynamical possibilities in multi-field inflationary 
scenarios are considerably richer than in single-field mod- 
els. An example is assisted inflation [l|, where a col- 
lection of fields evolve together, perhaps driving inflation 
that would not be possible if only one field were present. 
Configurations with multiple scalar fields have been con- 
sidered in the case of exponential potentials with equal or 
differing slopes. Later extensions of these models include 
the effect of curvature and a barotropic perfect fluid 3j , 
Bianchi backgrounds and polynomial potentials 5[. 
The stability of solutions as critical points in phase space 
was studied for decoupled [1, H, @] and coupled [1, [T(| 
exponential potentials, and for decoupled inverse power- 
law potentials [7J. Late-time dark-energy scenarios can 
be found in Refs. @, El- Multi-field scenarios naturally 
emerge in Kaluza-Klein compactifications Q and in M- 
theory, via the nonperturbative dynamics of a stack of 
M5-branes distributed along a S 1 /Z 2 orbifold [Hj]. In 
the latter case, the scalar fields and their exponential 
potentials are related, respectively, to the distances and 
instanton interactions between the branes. A model mo- 
tivated by strin g th eory is A-flation, where the scalar 
fields are axions [13j . 

Thus far, general conditions for the existence of clas- 
sically stable assisted inflation solutions have not been 
given in the literature. 1 This is in contrast to the related 
case of tracking quintessence models, featuring a single 
field and a barotropic fluid, where the tracking condition 
was found to be If 511 
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where V((f>) is the scalar field potential, 7 the fluid 
barotropic index defined by p = (7 — l)p, and prime 



a derivative with respect to the field value </>. The con- 
dition is often simplified to T > 1, as otherwise the field 
energy density does not grow relative to the fluid [l5| . 

The aim of this article is to derive those stability condi- 
tions. In the first instance we are interested in the case of 
n scalar fields <f>i each with an identical potential V{4>i). 
In that case Kim et al. [7| showed that, provided the 
fields all have equal value, the system was dynamically 
equivalent to a single-field model with field \ and poten- 
tial W(x) = nV(x/y/n)- This generalizes the correspon- 
dence for exponential potentials derived in Ref. [fj], and 
also shows there is no assisted inflation phenomenon for 
inverse power-law potentials |7j. However, while such so- 
lutions obviously exist, that paper did not address under 
what circumstances they are stable, which we do here. 
We then also consider some generalizations of this ba- 
sic set-up. None of our results require the Universe to 
accelerate, so they are valid also beyond the slow-roll ap- 
proximation. 



II. BACKGROUND EVOLUTION 

The background equations of motion with an Einstein- 
Hilbert gravitational action and n Klein-Gordon homo- 
geneous scalars read 



H 2 = — 
3 



W 



^ n—1 



3H(j) l + d l W = 0. 



(2) 
(3) 



where H is the Hubble parameter, k is the gravitational 
coupling constant, dots represent derivatives with re- 
spect to synchronous time t, W = W((f>o, <f>\, . . . , 4> n -i) 
is the potential with both self-interaction and interaction 
terms, and di — d/dtpi. For simplicity, we set k 2 = 1. 
The two equations can be combined to give 



H = - 



(4) 
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1 Inhomogeneous perturbations of multi-field configurations were 
considered in Ref. [13 . 



We are going to discuss the classical stability (that is, 
against homogeneous perturbations) of solutions of the 
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form 



i = 0, 



. , n - 1 , 



(5) 



i.e. where the fields all evolve together. This class of 
solutions includes multi-field fixed point solutions, e.g. as 
in Refs. 0, 0], but is much more general — the solutions 
need not be of scaling type (commonly defined as the 
scalar field kinetic and potential energies maintaining a 
fixed ratio, e.g. Ref. [ly|), and indeed need not even be 
inflationary. 

We begin by determining what forms of potential sup- 
port such solutions. Denoting derivatives with respect to 
</> with primes, this solution is possible only if 



diW\,,=V'(cl>) Vi,i. 



(6) 



where in the left-hand side <j)j — 4> for all j = 0, . . . n — 
1, and V is defined by this equation. If the fields are 
mutually decoupled, then each can be written with its 
own potential Vi obeying 



Vi(x)+Ai = V(x) 



(7) 



for all i. Here the Aj are constants, which can all be 
absorbed into some A = . Aj acting as a cosmological 
constant, so that W — J2i Vi{<fii) + A. We will see that 
this constant does not affect whether the solution is sta- 
ble or not, though it can alter the qualitative behaviour 
in the vicinity of the solution. Ordinarily we will consider 
A to be zero or negligibly small, so that all the fields have 
the same potential. 

If one allows for a coupling between the fields, it must 
still satisfy Eq. ([6]). This happens, for instance, when the 
interaction term is symmetric in all the fields, so that one 
can write the potential as 



W = J2Vi(<Pi) + f 



(8) 



for some functions / and g. Also, there can be more gen- 
eral situations, as the following three-field system shows: 
W = 0g + a(2b - a)- 1 ^- b cp b 1 + b{2b - a)' 1 ^. 



III. STABILITY 

Let us first consider the simple case of an expanding 
universe (H > 0) and no interaction terms, Vi — V. We 
do so using two different formalisms, the latter of which 
will prove the more powerful. 



A. Synchronous time formalism 

Perturbing the above equations of motion around the 
solution Eq. |5|), 

H(t) -f fT(t) + SH(t) & (t) -> cf>(t) + 5<f>i (i) , (9) 



and defining the (2n+l)-vector (T denotes transposition) 
X^iH^i^if, (10) 



the linearized equations of motion can be written in a 
matrix form SX = M5X, where the entries = (M)ij 
are 
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(12) 



The characteristic equation det(M — A/2n+i) = deter- 
mines the eigenvalues A = and 



~ 3H 



1±\ 1 



4(V" - 3n0 2 ) 
9IP 



(13) 
, (14) 



where A ± is n — 1 times degenerate. 

The solution is stable provided no eigenvalue has a pos- 
itive real part. Equation (fT3)) imposes the simple stability 
condition 



V" > , 



(15) 



while the real part of A~ is negative definite only if a 
stronger condition is valid: 



V" > 3n0 2 . 



(16) 



In the slow-roll approximation 3H 2 « nV, —3H(/> w V', 
this is nothing but the tracking condition 
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(17) 



but note that Eq. (| L6[) is valid also outside inflation. 

This looks like a natural conclusion exactly matching 
the expectation inherited from the known quintessence 
case. However in fact the interpretation is quite different, 
as can be seen by considering the implications in the 
simple situation of the power-law potential 



V 



r = i 
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(18) 



Clearly, this potential with p > (in particular, p = 2) 
does not satisfy the condition (fTB)) . implying Eq. (O is 
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not a slow-roll attractor, contradicting the known result 
fl7l ] that solutions for this potential do approach radial 
trajectories which, by field rotation, can be described as 
all the fields rolling together (see also Ref. [H[). 

The reason for this discrepancy is that this formalism 
fails to recognize that trajectories distinguished only by 
a shift in time, t — > t + St, are physically equivalent and 
should be identified, even if the apparent perturbation at 
fixed time appears to grow. This differs from the case of 
a tracking quintessence model, where the existence of a 
barotropic fluid acts as a "clock" , with reference to which 
time shifts in the scalar field evolution can be physically 
distinguished. 

The time-shift mode can be removed by defining quan- 
tities which are gauge-invariant with respect to the gauge 
group of gravitation (group of infinitesimal coordinate 
transformations) In the particular case of multi- 

ple scalar fields and homogeneous perturbations, this is 
8(j>i = S(j>i(t) — (j)St |l4[. Rather than writing down gauge- 
ready equations using the synchronous time formalism, 
we now employ a simpler alternative. 



B. Hamilton— J acobi formalism 

In the Hamilton-Jacobi (HJ) formalism 0, [2l|, the 
role of the clock is invested in a matter field, so that 
one can get rid of the unphysical degree of freedom 
represented by shifts in the time coordinate t. When 
many matter fields are present, one can choose one 
and study homogeneous perturbations on constant time 
slices. More specifically, we can take <fio — <fi as the new 
time coordinate and perturb only the other n — 1 fields 
4>i(4>)- The equation of motion for <f> is a constraint deter- 
mining the relation between different time coordinates. 
The first HJ equation is given by Eq. (@]), 



2H'{ct>) 



where 



(19) 



(20) 



here and in what follows i runs from 1 to n — 1 . 

(3H(p + W^d^, one gets n — 1 equa- 



Since d 2 = </> 2 <9? 



tions of the form 



W'4>\ + dA¥ = . 



(21) 



The expansion of the universe is encoded in the friction 
term and in <fi 2 via Eq. (fT!))) . Plugging Eq. (TT§|) in the 
Friedmann equation one gets 



2H' 2 - 3AH 2 + AW = . 



(22) 



When perturbing Eqs. (|2"Tj) and (I22p around Eq. ([5]), some 
simplifications occur, since <j)'l — and A = n on the 



background (8 A — 2 ^ S^'A. The final result is 

4> 2 84>'[ - V'Sti + Y^{d 3 d l W - djWy&j = , (23) 
3 

-H'SH' - 6HSH + J2 (v'Hi - 2 <50-) = , (24) 



where all background quantities are functions of <p. In 
general, the matrix djdiW is symmetric with all different 
entries. 

When the fields are decoupled, djdiW((f>) — V"((/))5ij, 



djW = 0, and j = i in Eq. (|2"3"|) . For simplicity we 
shall concentrate on this case. Again, the system can 
be written as SX' — MSX, where X is now a (2n — 1)- 
vector and the coefficients of the matrix M are given by 
Eqs. (J23J) and ([24]). Define the objects 



(25) 




(26) 



The eigenvalues of the decoupled system are Ao and the 
(n — l)-times degenerate pair X±. Concerning the first 
eigenvalue, it is normally convenient to take <j) > (ar- 
ranged if necessary by sending (f> — > — <f>) so that increas- 
ing <p corresponds to increasing time. The first eigen- 
value is then negative definite, as required for stability. 
If <j) were negative, then the two time choices would flow 
in opposite directions and the stability condition would 
become a positive-definite eigenvalue, the physics being 
unchanged. Therefore Eq. (f2"5)) holds in either case. 

With the convention (f> > one has that V < 0, as 
the situations of physical interest always have the field 
rolling down the potential and hence (j) should have the 
opposite sign to V. If V" > [Eq. (TTS])] holds there are 
two decaying modes related to A± (with the equality, one 
direction is flat) and Eq. |(SJ) is a classical attractor. If 
V" > V' 2 /4<^ 2 , the perturbations periodically overshoot 
the attractor. 

Notice that 



A 1 A_ = A1A* 



(27) 



meaning that the perturbations of the i-th field (i ^ 0) 
are related in both formalisms. However, they are not 
equal except in the slow-roll regime. To show this, let us 
denote as 5(f>^ the perturbation of the i-th field in the 
^-formalism. Using the diagonalized linear equations and 
Eq. (|27l), one finds 



The right-hand side is nonvanishing because 



(28) 



is not 



gauge invariant. However, it is negligible in the slow-roll 
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approximation. This is expected since, in this case, the 
eigenvalues would be almost constant and perturbations 
would be exponential: 5<fii ~ e x±t in the t- formalism, 
while 8(j>i ~ e A± ^ in the HJ formalism. Then the eigen- 
values would be individually related, A?j_ « <p\±. 



C. Interpretation 

The conclusion from these calculations is that the sta- 
bility condition for multi-field tracking in identical po- 
tentials is simply V" > 0, rather than a condition similar 
to the quintessence tracking condition. The interpreta- 
tion is in fact straightforward, and most easily seen by 
studying Eq. (|23[) for the perturbed scalar field evolu- 
tion. This equation is independent of the perturbed met- 
ric SH. Accordingly, a particular field perturbation is 
insensitive, at linear order, to the perturbation in the 
metric caused by any other field perturbation. The fields 
evolve independently of each other. Whether or not the 
field perturbation dies away, i.e. the fields approach one 
another, simply depends on whether those fields which 
happen to be further down the potential are evolving 
more slowly, which is the case provided the potential is 
less steep there. Assisted inflation does take place, but 
not because the fields are actually aware of each other's 
value and are drawn towards one another, but rather just 
as a generic property of trajectories on convex potentials. 

Multi-field systems which do not satisfy Eq. (|15P will 
never attain the configuration Eq. (JSj). On the other 
hand, Eq. (fT5|) guarantees the stability of the equal- 
field solution Eq. ([5|) to small perturbations, but does 
not tell us whether such a solution acts as a global at- 
tractor for the system. In the simple case of shallow 
exponential potentials, this global attractor property is 
known to hold [3]. However, it is easy to find counterex- 
amples. Consider the case of steep exponential poten- 
tials Vi(4>i) oc exp(A^i) for A 2 > 6. These potentials 
are so steep that there are no scaling solutions, as the 
potential energy falls off faster than the kinetic energy. 
Our particular solution is stable to small perturbations 
(V" = X 2 V > 0), but if the fields are initially set well 
apart the rapid decay of their velocity will prevent them 
asymptoting to equal field values. 

Note also that our condition V" > is a local stability 
condition referring to the position on the potential that 
the fields happen to be at during a given epoch. Some po- 
tentials may have V" with the same sign everywhere, and 
the stability or otherwise is then a global property of the 
potential. But potentials where V" changes sign may ex- 
perience sequences of epochs where the fields alternately 
converge or diverge from each other, the ultimate stabil- 
ity being determined by the sign of V" after its final sign 
change. 

We mention briefly the effect of lifting the assumption 
A = 0, which changes the value of H corresponding to a 
given location of the scalar fields (we assume we stay on 
the expanding branch). Looking at the eigenvalue equa- 



tion Eq. (fT3f . we see that changing the value of H can 
modify (or create) an imaginary part to the eigenvalues 
but can change the real part only by an overall multi- 
plier, leaving its sign intact. Hence introducing A can 
change the type of attractor (oscillatory or monotonic), 
but not whether the potential gives stable solutions or 
not. A positive A can also impact on whether or not 
the solution Eq. §| is a global attractor, in the sense 
specified above. Such cosmologies will be asymptotically 
de Sitter, with rapidly-decaying field velocities that may 
prevent approach to the equal-field solution. 

The impact of A is somewhat less transparent in the 
Hamilton-Jacobi formalism, Eq. (|26|) . where at first sight 
A± seems independent of H even though we have verified 
that within slow-roll A!j_ « cj)\±. The resolution is that 
at a given location on the potential, the introduction of 
A will modify the value of <fi in Eq. ([26]) . again poten- 
tially introducing or modifying an imaginary part to the 
eigenvalue. 

One might ask how our conclusions would change if 
a fluid component with equation of state p = (7 — l)p 
were added, to give a multi-field quintessence scenario. 
Adding a new component changes the nature of the dy- 
namical system, and the new degree of freedom may well 
induce an instability, even if the fluid component is taken 
as initially sub-dominant. E.g., for sufficiently-steep ex- 
ponential potentials the late-time global attractor has the 
field scaling with the fluid Nevertheless, one would 
expect that within the field sector the stability condition 
remains the same, i.e. provided V" > 0, the stable solu- 
tion will still maintain equal field values. It is easy to see 
that this is true if p = in the background solution, pro- 
vided 7 > (the extra degree of freedom would decouple 
from the scalar fields with eigenvalue oc —jH/(p). 

Notice also that when 7 = 2 (stiff matter, p = p) 
Eq. |T]) reduces to the pseudo-tracking condition Eq. ([±7]) 
found in synchronous time for slowly-rolling fields. Stiff 
matter decays as p ~ a~ 6 , this being the fastest rate 
at which a scalar field density can be diluted (kinetic 
regime, <j) 2 3> V). Since slow rolling prevents this con- 
dition, for all purposes the fluid is negligible relative to 
the scalars, which is precisely the configuration leading to 
Eqs. (fT6|) and IT7|) . However, this implies that one cannot 
use the barotropic fluid as a reliable clock throughout the 
whole evolution of the system, and the calculation in syn- 
chronous time becomes nonphysical. In other words, clas- 
sical stability in multi-field cosmologies with a barotropic 
fluid can be consistently studied in synchronous time for 
7 < 2 but not in the limit 7^2 (no fluid), where the 
HJ formalism is more appropriate. 

IV. GENERALIZATIONS 

A. With field couplings 

The next simplest case is that of a symmetric interac- 
tion, Eq. ([5]), so that the nondiagonal entries of djdiW 



5 



are all equal, and we can define the useful quantity 

y{ff>) = d j d i W = d j d i f ) j^i. (29) 

Defining 



v 

20 2 



1±A 1 



, V" + (n- 2)y 



yr. 



(30) 



one can show that the HJ eigenvalues are Ao and the 
(ri — l)-times degenerate pair Stability is achieved 

if 

V">y. (31) 

We can also consider the case where the coupling is sym- 
metric in all the fields except one which remains uncou- 
pled (djW = 0). It is natural to choose the decoupled 
field as the clock field <fi — <fro (though one would be free 

to choose another). Then the eigenvalues are Aq, A 



O) 



(i) 



and the (n — 2) -times degenerate pair A ± 
bility condition is 

V" > max{y, -(n-2)y) . 



while the sta- 



(32) 



With interactions included, stability now requires 
steeper self-interaction potentials, which is a confirma- 
tion of the claim made in Refs. [H, Q that multi-field 
inflation with cross-couplings is less likely to happen. If 
y < 0, the larger the number of fields the steeper the 
coupling term. As an application of Eq. (|32p . consider n 
scalars with potentials Vo = Ccj) p , Vi — (j) p , and power- 
law interaction / = B{(f>i ■ ■ ■ <^ n _i) p /( n_1 ), where B, p 
and C = 1 + B/(n — 1) are constants. Let also B,p > 0. 
Then y = Bp 2 <j}P~ 2 / (n — l) 2 and the solution is stable if 
p>C. 

If the cross-coupling is asymmetric, the symmetric ma- 
trix djdi W((f>) and the eigenvalues of the linearized equa- 
tion are more complicated; however, for a given model 
it is straightforward to check its classical stability via 
Eqs. (E3J) and (EJ). 



B. More general backgrounds and the braneworld 

We can generalize all the above results to other back- 
ground equations and solutions. It is instructive to con- 
sider background solutions of the form 



(33) 



for some real cti and ft. Equation (|33|) includes the solu- 
tion Eq. (fTI))) (a, = 1 = ft), scaling solutions [2, [3, the 
power-law inflationary attractor studied in Ref. |7|, and 
other situations where, for instance, some fields are sub- 
dominant with respect to others (ai rs for some i). The 
condition on the potentials for Eq. (|33|) to be a solution 
is 



d,W\ 



V'W + il-pi)- 



(34) 



Again, V' is defined to be the potential of <po, differen- 
tiated with respect to (fio, where all fields are set to be 
equal to 4>o ; in other words, Eq. ([6]) with i = 0. This 
condition is fulfilled, e.g., by slowly-rolling fields with 
different polynomial potentials. On this background, 
A(<j>) = 1 + EjKft ) in Eq. CED). 

In high-energy braneworld scenarios inflation can be 
sustained by steeper potentials [13] and it would be in- 
teresting to see how the stability conditions are affected. 
Neglecting brane-bulk energy exchange and the contri- 
bution of the Weyl tensor, the Klein-Gordon equation 
© [as well as Eq. (fM)) ] is unchanged, while the effective 
Friedmann equation on the brane at high energies reads 



(35) 



Here, 9 = in the pure 4D regime, 9 = 1 in the high- 
energy Randall-Sundrum braneworld, and 9 = —1 in the 
high-energy Gauss-Bonnet scenario; Kg is the effective 
gravitational coupling felt by the matter (with energy 
density p) on the brane. See Refs. [23| for a review of 
this formalism and references to the original results. 

We continue to work with the more general background 
solution given by Eq. (|3"3"|) . In units k 2 , = 1 — 9/2, the HJ 
equations become 



H' z = 



2H' 

~Wa 

AH 26 



6H 2 



W 



(36) 
(37) 



Linearizing Eq. (|37|) one gets 
SH' = 



^(3 + 6»e) SH 



+— r- (2 diWSfa - 4> 2 5A 



with 



and 



H 



Acf) 2 

e = ~~H 2 ~ m 2 ^ 1 



5A = 2j2<*ifc 



(38) 



(39) 



(40) 



Also, since </)" = a, ft (ft - l)0 ft " 2 ± 0, the perturbed 
Eq. (|21[) now contains the metric perturbation: 



12H 1 



Ac)) 2 



V 



A 



(41) 
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The eigenvalue equations can be solved numerically or 
analytically for particular cases. A simple possibility is 
Pi = 1 for all i, whereupon <ft" — and many of the terms 
in Eq. (|41[) vanish. The background condition on the to- 
tal potential becomes diW — ctiV', and the perturbation 
of the Hubble parameter decouples from the scalar ones. 
In this case, the form of the effective Friedmann equa- 
tion affects only Eq. ([55)) . The eigenvalue Ao is positive 
definite only if [see Eq. (|38]>] 

- Oe < 3 , (42) 

which is automatically satisfied for 8 > (4D and 
Randall-Sundrum braneworld). In the Gauss-Bonnet 
scenario, it must be e < 3, which is true during accel- 
erating expansion. 

When there is no coupling between fields (djdiW = 0), 
the other stability condition Eq. (fTS"|) becomes 

onV" > . (43) 

For quadratic potentials this states that none of the 
scalars is tachyonic. In the presence of interactions and 
for a>i ^ a, , the degeneracy of the eigenvalues is broken 
and more involved relations arise. 



tions. These can be summarized as follows. 

• For the solution Eq. ([5]), if the fields are nonin- 
teracting the only condition is V" > 0, while for 
symmetric interactions with positive coupling the 
self- interacting potentials must be steeper. 

• For the solution Eq. (|5o| with Pi — 1 and no cross- 
coupling, the stability condition reads aiV" > for 
all i. 

• These results are valid for both four-dimensional 
cosmology and high-energy Randall-Sundrum 
braneworld. On a Gauss-Bonnet braneworld, the 
extra condition e < 3 is required. 

• In all other cases, the linearized dynamical system 
can be studied numerically via Eqs. (|38|) and (|4Tj) . 
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